arXiv:1506.09097vl [physics.atom-ph] 30Jun2015 


Control of atomic transition rates via laser light shaping 
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A modular systematic analysis of the feasibility of modifying atomic transition rates by tailoring 
the electromagnetic field of an external coherent light source is presented. The formalism consid¬ 
ers both the center of mass and internal degrees of freedom of the atom, and all properties of the 
field: frequency, angular spectrum, and polarization. General features of recoil effects for internal 
forbidden transitions are discussed. A comparative analysis of different structured light sources is 
explicitly worked out. It includes spherical waves, Gaussian beams, Laguerre-Gaussian beams, and 
propagation invariant beams with closed analytical expressions. It is shown that increments in the 
order of magnitude of the transition rates for Gaussian and Laguerre-Gaussian beams, with respect 
to those obtained in the paraxial limit, requires waists of the order of the wavelength, while propa¬ 
gation invariant modes may considerably enhance transition rates under more favorable conditions. 
For transitions that can be naturally described as modifications of the atomic angular momen¬ 
tum, this enhancement is maximal (within propagation invariant beams) for Bessel modes, Mathieu 
modes can be used to entangle the internal and center of mass involved states, and Weber beams 
suppress this kind of transitions unless they have a significant component of odd modes. However, 
if a recoil effect of the transition with an adequate symmetry is allowed, the global transition rate 
(center of mass and internal motion) can also be enhanced using Weber modes. The global analysis 
presented reinforces the idea that a better control of the transitions between internal atomic states 
requires both a proper control of the available states of the atomic center of mass, and shaping of 
the background electromagnetic field. 
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I. INTRODUCTION. 

Modifications of the internal states of an atom due to its coupling to the electromagnetic (EM) field are, under 
standard conditions, highly restricted. Once the proper selection rules are satisfied, the enhancement or inhibition 
of transition rates can be based on the control of the evolution of the atomic center of mass, and the appropriate 
selection of the features of the surrounding electromagnetic field. The latter includes its frequency spectrum, angular 
spectrum and polarization. In this work, we make a systematic study of the feasibility of modifying atomic transition 
rates by tailoring the electromagnetic field, given different schemes for the center of mass dynamics. 

In general, previous analysis that had the enhancement or inhibition of internal state transition rates in mind, can 
be classified as those which study atoms confined in electromagnetic cavities [I] or nearby surfaces (2j, and studies 
in which laser light is shaped in open space to create electromagnetic fields with a predetermined structure m 3 . 
The first scenario, in general, gives rise to polarization restrictions: the natural modes would be, e. g ., transverse 
electric or transverse magnetic EM fields satisfying the adequate boundary conditions on the cavity surface. In the 
last scenario, most descriptions usually assume circular or linear polarization with respect to the main propagation 
axis of a given laser beam. Here we present a general analysis that includes both cases. 

We take into account the relevance of the central field description of electrons in atoms, and the fact that EM fields 
can be described as superpositions of modes with well defined total (orbital plus polarization) angular momentum. 
The accompanying symmetries of both atom and EM field are essential to get a proper description of the transition 
rates under realistic conditions. In all the calculations, the role of the atomic center of mass will be emphasized. 

In Section II, a general accurate description of the so called scalar and vectorial EM modes via their wave vector 
spectra is presented; using it guarantees the validity of our forthcoming analysis beyond the limitations of the usual 
paraxial approximation. In Section III, the transition matrix elements of the minimal coupling interaction Hamiltonian 
of a nonrelativistic hydrogen like atom are written in terms of products of a center of mass factor, an internal atomic 
state factor and an electromagnetic factor which couples the former two terms. Based on this, a modular description 
of transition rates involving all the degrees of freedom of the system atom-EM field is developed. In Section IV, we 
revisit the general properties of the internal transition factor and the standard classification in terms of the coupling of 
different atomic multipoles to the EM field. We also show that the developed formalism can be used to study internal 
transitions which require a relativistic description without modifications on the center of mass and electromagnetic 
factors. In Section V we stress the main features of the center of mass term when it is described by a quantum wave 
function. The limitations of such an approach under realistic conditions are also briefly discussed. In order to compare 
with standard studies of atomic recoil effects, in Section VI, we apply our formalism when the center of mass states 
is described in terms of plane waves. In Section VII, we report the explicit calculation of the electromagnetic factor 
for spherical modes, Gaussian modes and propagation invariant modes for a center of mass motion with spherical 
symmetry. This could correspond, for instance, to atoms in free space or in a spherical trap. In section VIII, an 
optimization procedure to enhance transition rates by selecting the proper parameters of a Laguerre-Gaussian beam 
or of a Bessel beam is explicitly worked out; the trap that confines the cold atoms is assumed to have circular cylinder 
symmetry. Finally, some conclusions derived from this study are given in Section IX. 


II. THE LIGHT FIELD: STRUCTURED BEAMS AND THEIR ANGULAR SPECTRA. 

The electric field E(r, t) of an electromagnetic wave in otherwise free space can always be written as a superposition 
of plane waves with wave vector k, frequency w and polarization b: 

duSioj) f d 3 ke^-^a^l - ^ 2 /c 2 -fci)e^(k x ,k y ;u;). (1) 

b K 

In this equation S(uj) denotes the frequency spectra, the delta function guarantees the fulfillment of the dispersion 
relation, k± = ^Jk% + ky , and d b ^(k x ,k y ',u>) is the vectorial spectrum in k space for the polarization 6; it is charac¬ 
terized by a set of parameters n that are required to specify each electromagnetic mode of frequency w. The magnetic 
field associated to the EM wave is given by 

bm = EE/ dboS(uj) f d 3 ke l( - kr ~^6(\k z \ - ^uVc J -ii)i8<‘)(kx,k y ;w) (2) 

b K ** ** 


with 


<8(. b) (fc x , ky\uj) = (c/w)kx £^\h x ,ky\u)). 


( 3 ) 
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In the Coulomb gauge, the electromagnetic vector potential for each mode is 

Af>(r ,t) = j d 3 ke^ k ' r_w ^A^(k; w) (4) 

= d 3 ke^ r ~^S(\k z \- y/u>*/c*-k a ± )<£W(k x ,k y] u>). (5) 

The polarization b reflects the transverse structure of the electric field V • E(r, t) = 0, so that, 

k • <£W = 0. (6) 


One selection of the resulting two independent polarizations is defined in terms of a preselected vector a as follows 
(from now on, v denotes the unitary operator along v): 


£l E (k x ,k y -,uj) = £ EE [k x a }¥ EE \k x ,k y ), 
&l M {k x , k y ; w) = £ J M [k(k ■ a) - a}¥™ ] (k x , k y ). 


( 7 ) 

( 8 ) 


These modes are the well known transverse electric (TE) and transverse magnetic ( TM ) modes frequently used in 
cavity electrodynamics (although they can also be generated in free space [8]). The vector a in such systems is usually 
chosen as an axis of symmetry of the cavity. The normalized functions V I / [vL (k x , k y ) are known in the literature as the 
angular spectra of the mode, the factor £W is proportional to the mode amplitude. Specific examples of (t'S^kx, k y ) 
are given below. If a is parallel to the z axis, the beams are said to be vectorial EM structured beams with a privileged 
propagation axis along the z-axis; with this selection 


a — e~ = eo 


and e± = (e^ ± ie y )/V 2, 


( 9 ) 


and 


e 

,(tm) = 


{TE) Ekxe. = -^ £ r(e* Vk e_ - p"^ k 


= k(e ; ■ k) - e z = 


y/2cui 
k i k z c' 


y/2cc 


(e" 
-(e* Vk e_ 


e+) 


k 2 , c 2 . 


00* 


( 10 ) 


Here is the azimuthal angle in the wave vector space, e ±lLpk = (k x ± ik y )/k±. 

Another selection of polarization basis originates from the paraxial description of light beams. In this case, the so 
called scalar structured beams are described as follows. First, a fixed direction is selected a 3 , and two other orthogonal 
directions are identified a. 2 , and ai. Then, 


<B 1 K (k 1 ,k 2 ;uj) =£*[ai - (fei, k 2 ) 

k • a 3 

£% u (ki,k 2 ;u)) = £% [a.2 — ^4^a 3 ]§^(fci, k 2 ) 

k • a 3 


£le^l(kiMl 


( 11 ) 

( 12 ) 


In the paraxial limit a 3 usually corresponds to the main direction of propagation of the beam. To be specific, let 
it be the z-axis, i. e., a 3 = e z . The paraxial condition \k z \ ~ oj/c yields (k • a.j)/(k • a 3 ) < 1, i =1, 2, so that the 
component of the field along a 3 is usually neglected; if a 3 = e x and a 2 = e y , a 3 and a 2 define the two orthogonal 
linear polarization directions. 

Scalar EM modes may, under adequate circumstances, be approximately generated using a monochromatic laser 
beam with the proper polarization that impinges a space light modulator or a structured diffraction plate. To that 
end, each optical device is respectively programmed or designed to guarantee the appropriate boundary condition 
that defines each mode at the plane z = 0: 


a, • E, 


t,u(x,y,0) = E/ d 3 ke < <*- s+ *»i'>J(|k,| - yju?/# - fci)a, • <EW(k*,k y . u ) 


d2 ke < ( fc - a + fc »»)5‘$W h) (k a ,k I ,) 


^ £^%{k x ,k y ) = jdxjdy e - i ( fc - x+fc »*')a i ■ V KiU (x,y, 0). 


In Eq. (13), C w / c denotes the interior surface of a circumference of radius lo/c in the k^ space. 


(13) 

(14) 
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Vectorial EM modes can be written (an generated) as superpositions of scalar EM modes and viceversa by the 
selection of the proper angular spectra; taking a = a 3 , 


e ( TE ) = k x a 3 = k 2 [&i ~ (h/h)& 3 ] ~ &i[a 2 - = & 2 e (1) - Aqe (2) , 

e (TM) = k(a' 3 • k) - a 3 = &ife 3 [ai - (fci/fc 3 )a 3 ] + fc 2 fc 3 [a 2 - (/c 2 //c 3 )a 3 ] = fcifc 3 e (1) + fc 2 A; 3 e (2) . 


(15) 


Notice that the superpositions depend on k. Quantization of structured modes in terms of vectorial beams [9] is 
usually more direct than with scalar modes due to the evident orthogonality between TE and TM modes, in contrast 
to the overlap between the scalar modes defined by Eqs. fllipl i- 

The purpose of this work is to present a general formalism to compare, in a quantitative way, the dependence of 
atomic transition amplitudes for different selections of polarization and angular spectra of a coherent light source. 
For the selection a = e z = a 3 , explicit results will be given for the following angular spectra : 


• Gauss modes with waists w x and w y : 


mm;J k X ,ky) = 


—w tkt —vj ,, Ar 


(16) 


• Laguerre Gaussian modes with waist w Q , topological charge m and order p : 

p 


E(-!) a ( P p + - 1) (1?) 

a —0 ' ' 


• Propagation invariant beams with analytical expressions. These beams, under idealized conditions, maintain its 
transverse structure with respect to a given direction that here we take as e 2 . They have the angular spectra 
structure 


= A/-W 


k 1 


(18) 


Explicitly: 

— In circular cylinder coordinates, they are known as Bessel beams nun]: 


\T/( b ) (h b ) — \f( b ) K±) 


(19) 


— In elliptic coordinates with an ellipticity parameter q and order m, they are known as Mathieu beams 

mm-- 


,m;a) ^ ( Ce m(^ g) + ^e TO (y? k ; q)). 

where ce m (ip; q) and se m (tp; q) are Mathieu functions of order m in standard notation |13j . 

In parabolic coordinates, they are known as Weber beams mm- For even parity and aperture a: 


( 20 ) 




S(k _L - K_l) e i«ln|tan V k/2| 


k _L 


V 


7T| sin (fk | 


( 21 ) 


Odd parity modes are obtained by replacing | sin in the denominator by —i \ sin<^k|- 

Each of these elementary modes can be combined to yield, for instance, even standing transverse Bessel modes 
l,m;u(k x , k v ) + -m-u(k x , k y ) or traveling Weber modes by superposing even and odd modes with a ±7r/2 phase 

factor. The experimental generation of approximate propagation invariant beams can be implemented by replacing 
the delta factor by a Gaussian function of k±, centered around and with a width much smaller than 

The calculation of the normalization factors Wifi, is discussed in Appendix A. Their selection is done to allow an 
equitable comparison between the effects of different modes in the atomic transition rates. 

We shall also study the results obtained for vectorial spherical waves. In this case, the vector a that defines the 
polarizations is chosen as the gradient operator in k space, Vk and it is applied to the spherical harmonics in that 
space. Spherical electromagnetic waves are described in detail in Appendix B. 
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III. ATOM-LIGHT INTERACTION. 

A. Coupling of the electron current to the electromagnetic field: the Hn interaction Hamiltonian. 


For simplicity we consider a hydrogen like atom in conditions that do not require a relativistic treatment. In the 
next section, we describe how these restrictions can be relaxed to include multi electron atoms and relativistic effects. 

We are interested in the evaluation of the transition probability by which an ultra cold atom, initially described 
by a wave function that can be written as a superposition of functions that are products of a center of mass function 
$ k cm(R) and an internal function </> K rei (r re ;) with quantum numbers k cm and n rel respectively, 

(R)0 K rei (l re ;), 

changes its state to a wave function with the same structure, that is a superposition of wave functions of the form 

'rei) = $ K CM(R)^ K „i(l> e ;), 

and simultaneously absorbs a photon described by a mode of the EM field A^\ with polarization b and structure 
labeled by k. In the former expressions the position of the electron is given by the coordinate 

r e = R- + (22) 

where R refers to the atom’s center of mass, M e is the electron mass, fi = M e M^ / (M e + Mm) is the reduced mass, 
Mm is the nucleus mass, and r re i is the relative coordinate. <I> k cm(R) and </> K r-ei (r re i) are eigenfunctions of the 
corresponding Schrodinger equation with initial eigenvalues Eq^ m and E^j , and final eigenvalues E^ 2r and E^. 

We focus - mostly- in the case that the transition is induced by the interaction Hamiltonian 


Hi i = -JfPrei ■ Af(r e 


(23) 


with q e , the electron charge and p re i the relative momentum. This Hamiltonian can be understood as a result 
of the coupling of the electron current q e p e /M e to the electromagnetic field A^(r e ) in the approximation where 
p e = ( M e /(M e + Mm ))P + p rei ~ p rei valid whenever Mm M e ( which is satisfied under standard conditions). 
Corrections to this approximation are described below. 

If A^ satisfies Coulomb gauge ( V ■ A ^ = 0) within first order perturbation theory, the transition amplitude T 0 f 
between the elementary wave functions ^(R, r re ;) and \I/f(R, ?rei) is given by 


Tqf = (*f\Hii\* 0 ) = Jr(£<rf ~ E< fJ) J d 3 r rel d 3 R[^* f { r rei , R)r re ^ 0 (r, ei , R)] 


^ e M p A " )(R+ M/ rei) ]’ 


(24) 


in the length form. Taking into account Eq. ([5]), 

A?>(r e ) ==™J d 3 ke i( - k < n+ ^ r ^-^S(\k z \ - ^/c 2 - fc 2 )(*W(fc x , k y] u). (25) 

The structure of Eqs. |7]|8l and Eqs. ( TTfl2 > let us write the absorption transition amplitudes Top in terms of the 
light angular spectrum as follows 

duS{u)e~ iAt J d 3 k S(k z - ^ - fci) Wk; i4 M )J rel (}r, «5 d ; K r /) ■ e^¥^(k x , k y )( 26) 


Tof = Vrei — 


Icm{ k; Kq 

J rel( k ’ 


GM. k CM\ _ 


,rel. ,,rel 
’'O , lx F 


F W 

>■/ 


(R )e ik,fl icM 


d' R&* cm 
d 3 r rel <j)* l (r rel )e 


(R) 


+'W k ' r '-»l 


(*Vez) 


TJrel — 


Qe^^rel 


rei — 




UCM = 


^CM ^• 


( 0 ) 

CM 


A = W — UJ re i — UJCM- 


(27) 

(28) 

(29) 


M e h h 

The corresponding emission amplitudes induced by the presence of the coherent EM field Tfo satisfy the relation 
Tfo = Tq F . Eq. (29) decomposes the transition amplitude Tof in factors that involve separately the internal degrees 
of freedom of the atom J re i, and the degrees of freedom of the center of mass Icm- These factors couple, in a modular 
form, to a coherent EM field via its angular spectra ifr$ui(k x , k y ), and its polarization e ^ (which in general depends 
on k in a way determined by n and b). 
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B. Other interaction terms resulting from the nonrelativistic minimal coupling Hamiltonian. 

The non relativistic Hamiltonian of two charged particles (the nucleus and the electron) interacting with the light 
electromagnetic field via the minimal coupling scheme, involves, besides the Hamiltonian Hu, two terms. One of 
them, for a monochromatic wave at resonance A = 0, gives rise to the transition amplitude 

(^ f \Hi 2 \^o) = itocM [ d 3 r rel d 3 R[** F (r rel , R)R^(r rei , R)] [q e A«*(R + ~^~r rel ) + q N Af (R-r rei )] 

J ivi N 

- J d 3 r rel d 3 R[<H* F {v reh H) r rei $;(r rel ,R)]A«*(R-i-r rel ). (30) 

with qN the nucleus charge and rjy = R — (/z/M/v)r re ; its position. The first term, in the last equation, results 
from the difference between the electron momentum and the momentum of the relative coordinate in the electron 
current-EM field coupling, the second and the third represent the coupling between the nucleus current qNPN / Mn 
to the electromagnetic field A^(rjv). 

For atoms trapped in a harmonic potential with characteristic frequency Cl, the center of mass motion can be 
naturally measured in terms of R ~ \JH/CIM F while the transition frequencies are integer multiples of u>cm ~ Cl, 
thus locmR ~ s/hCl/Mr which corresponds to ~ 0.1 m/seg for a harmonic trap with Cl ~ lOOKHz and a hydrogen 
atom. For the internal motion oj re i ~ Z 2 a 2 sc nc 2 /2h and r re i ~ h/Zp,caf sc , with af sc the fine structure constant; so 
that UJ re ir re i ~ Za/ sc c/2. Thus, 


q e uj re ir re i/q N uJcMR ~ a f sc c^/ M T / hCl ~ 10' > 1. 

That is, the natural units to measure the matrix element 

UrelQe (<j) K CM \ R \ <j> K CM ) (^rei | V re l \ <j> K re, ) , 

that defines the transition amplitude To.f, are much bigger than those related to 

UCMQN (^«g M I Re lk ' R | <j) K CM ) {(j) K rel | e ± * k ' rreI | <j) K rel ) . 

Since, under standard conditions /i <C M F and q n > q e , the order of magnitude of the relevant transition rates 
induced by any term in Hj 2 is expected to be much sm aller than its equivalent one induced by Hu. 

Note, however, that transition amplitudes of H 12 , Eq. (301, can also be written in a modular form similar to Eq. (291 
that involves the vector integral 


•CM 


_ik-R,T 


(k; K ^ m - k£ M ) = J d 3 R$* K cM (R) 

= -*V k / CA f(k;K^ M ;K^ M ) 


e“" *-R$ k cm(R), 

CM. ..CM\ 


and the scalar integrals 


J r ± (k;«5 rf ;«^) = J d 3 r re i<t>* K rei {r re i)e ±i ^ k ' rr,>l <j) K rei (r re ;), 


for the first two summands in Eq. (30), and the scalar integral Icm , Eq- (27), and the vector integral 


J re;( k b 


.rel. „rel 
"0 J ^F 


) = j d 3 r rel ^ rel )e (r re j) = J^( 


_U. ^rel. ,.reZ\ 
K, /Aq , 


= —*V k J rei (k; 




! ) 


(31) 

(32) 

(33) 

(34) 

(35) 


for the third summand in Eq. (301. So that the main features of the modular evaluation of the transition amplitudes 
developed in the next pages can also be applied to Hj 2 - 
The third term Hj 3 in the minimal coupling Hamiltonian 


H I3 = 


2 


\ A ( r e )\ 2 + ^-\ A ( t n )\ 2 


2 M, 


N 


(36) 


is usually interpreted in terms of two photon processes: introducing a complete set of intermediate center of mass 
and relative coordinate states the corresponding transition amplitudes are (g e ,Jv/2M ei jv) ^/(^fI (re.jv)!’^/) ■ 
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(’J'/IAiV, jv)|’I'o)- These transition amplitudes are proportional to ot 2 sc and, thus, generally smaller than the 

relevant ones arising from the Hj\ Hamiltonian. Note that the H 13 matrix elements can also be written in a modular 
form involving the scalar integrals J pel and Iqm defined above. 

It is important to mention that working with the minimal coupling scheme in the Coulomb gauge gives equivalent 
results [16. to those obtained in the PZW formalism C 21 118] . This formalism is obtained by applying a gauge 
transformation to a minimal coupling Lagrangian. It makes evident the relevance of the magnetic and electric dipole 
moments of the neutral atom evaluated with respect to the atomic center of mass; the electric and the magnetic fields 
of the photon field appear directly in the resulting interaction Hamiltonian which also includes a term quadratic in 
the EM field. Detailed studies of transition probabilities of atomic systems have been performed using the linear 
electric field part of the interaction Hamiltonian of the PZW formalism for a Laguerre-Gaussian beam in the paraxial 
approximation HIS 1 - 


IV. INTERNAL STATE FACTOR J r 


Notice that the matrix element J+,, Eq. (281, contains the information of the final and initial internal states of the 
atom and corresponds to that worked out in standard atomic physics calculations. From it the matrix element J~ ei , 
Eq. ( [34| , can be calculated. In this section we revisit key features of these matrix elements that are essential to follow 
the forthcoming calculations. 

Whenever the internal state of an atom is properly described by a central field interaction, that is, if tp Krd ( r ) = 
Rnim{r)Yi m ( r) (with Yj m the l, m spherical harmonic) the internal state quantum numbers are n re i = n, l, to, and a 
relevant relation that makes easier the evaluation of the relative coordinate integrals J re i is 


e * k r = E i l j^kr)Y^ m( (y)Yl mt {r), 

ft,me 


(37) 


which expresses plane waves in terms of spherical waves, in this equation je is the spherical Bessel function of order 
l. Each term in this series is an eigenfunction of the orbital angular momentum operators, L 2 = (r x V) • (r x V) and 
L z = (r x V) z . If Y) j?71( ,(k) were taken as an angular spectra and the proper polarization vectors were used, vectorial 
spherical electromagnetic modes which carry orbital angular momentum would be obtained as described in Appendix 
B. Using Eq. (371, it results 


J+,(k;4 e V/) = E / ro '’^(FM/e;C';4 ei )e;/i, m( (k), 


(38) 


tf,m e ,m ph (nu/cMe, K r 0 el ;K r /) = (^ w / cM e)w (4) (l, l', £, U m a , m F , 7%, m ph ), 

with e _ mph , vriph = 0, ±1 the polarization vectors defined in Eqs. (|9j, 

POO 

K rei K rei(^u/cM e ) = / drr 2+l R* nolomo {r)j e {fj,u}r/M e c)R nF i FmF {r), 

’ ” ° ’ F Jo 

and the angular integral 

t; (4) ( l 0 , l F , £, 1; m 0 , m F , me, m ph ) = 4? ri e j dtt r Y e * me (r) Y z * mo (r )Y lFmp ( r)Y lmph (r). 

The dependence of J)( ei (k; re™ 1 ; n r p l ) on k is through its modulus w/c and its orientation Ye t m t { k) as a manifestation 
of the Wigner-Eckart theorem. 

Particular useful forms of R n i m {r) correspond to linear superpositions of Slater type orbitals (STO) e _/3r r n_i ” 1 . 
For them, the relevant integral can be directly computed 


POO 

Qpa+0 F -,u, = J o dre-^ 0 + ^ )r r a je{^r/M e c) = 2t+1 ^ 


T(£- 


1) 


+ p F ) s+1 r(£ + 3/2) 


( yu Y / l+ 5 + 1 i + s + 2 / (UV 

\cM e (/3 0 + (3 F )J 2 2 ’ 2 ’ + 7 ’ \M e c(P 0 + /3 F )J )' 


(39) 









For a hydrogenic atom /? = (Z/n)(af sc nc 2 /he). Since the energy associated to changes in the center of mass state 
of an atom involve in general energies much less than those of electronic states, close to resonance w ~ uj re i ~ 
(. Z 2 a 2 sc ^c 2 /2h)(l/n 2 — l/n F ) and the term 


IJM ^ ZafscH ( 1 _ 1 _\ 

M e c(f3 0 + /3 f ) 2 M e \n 0 n F )' 

Note that, unless Zaf sc is not too small, i. e., for atoms requiring a relativistic description of their internal state, 
the series that define 2 F 1 can be cut at the zeroth order term 2 F 1 ~ 1. That is, under non relativistic conditions, 
retardation effects of the atom-light interaction can be neglected, and the internal atom transitions in terms of static 
electric and magnetic multipoles is valid. Besides, for Za.f sc -C 1, Q^ +p F - u rapidly decreases as £ increases: the 
highest possible values of 3 e ’ mt ’ rnph (/rw/cM e ; n r 0 el \n r p l ) can be found for £ = 0 , that is, electric dipole transitions 
induced by the Hamiltonian Hji . For such transitions, the companion spherical harmonic loo does not depend on 
k. The next order term in the series of J^ e; (k; k f ) depends on k through F l m (k). If the quantization axis of 
the atom coincides with the z-axis, 


n,o(k) 



*i,±i(k) 


f ~3~k±c ^ ±i(pk 

V 87 T U! 


(40) 


This term can be related to electric quadrupole transitions and magnetic dipole transitions induced by the Hamiltonian 
Hu. However the latter would not be properly described unless electron spin effects had been incorporated. 

The angular structure of given by reflects the vectorial character of the electromagnetic field through 

the presence of the Yf - m , iih factor, and the relevance of angular structure of that field through the orbital angular 
momentum term Y^ mi . Both of them can be coupled to yield an effective total angular momentum 


Y,- 


Yhlph 


(r)Y^ 


:( f )= Z! 

(-T ,rriT 


(2£ + l)(3)(2£ r + l) 


47T 


1 £ 

-TTlph, 'mt 


It 

rriT 


1 £ 
0 0 


Yu 


(?) 


(41) 


A similar expression can be used for Y* mo Yi F:mF so that the angular integrals v^ have a closed expression in terms 
of the 3 -j symbols. The selection rules for each £ multipole result from the structure of 3 -j symbols. 

For the evaluation of the transition amplitudes related to the Hamiltonians H 12 , Eq. (30), and H 13 , Eq. (36), it is 
necessary to calculate the internal scalar transition factors J Fel , Eq. (33). Using the series expansion given by Eq. (371, 
these factors can also be written in a completely analogous form to that of the vector factors J± ei . Now the functions 


and l F , £, 0; m a , mp, m #, 0) are the relevant ones. Whenever retardation effects can be neglected 

the standard condition kr re i << 1 makes that the term £ = 0 in the series expansion of the scalar factor J Fel gives 
a negligible value, and the terms with lowest £ value that may give a significant contribution correspond to £ = 1 , 
mg = 0, ±1 instead of the £ = 0 term of the electric dipole transitions induced by Hu. That is, the scalar integral J ± 
term associated to dipole transitions has the structure J± ipole rei ~ 93^i/ 4 )Yi im( ,(k). A similar argument follows for 
higher multipoles. 

For multiple electrons atoms, most of the equations given above are still valid both in the context of the independent 
particle model with a central field approximation, and for treatments that yield natural orbitals m after including 
electronic correlation effects. 


A relativistic description of the internal state evolution of a hydrogen like atom can be done in terms of the solution 
of the Dirac equation. Then the effective interaction Hamiltonian for the electron-EM field coupling could be taken 
as Hp = whenever the CM motion is not relativistic, with 7 ^ the Dirac matrices and the electromagnetic 

potential evaluated at the position of the electron undergoing the atomic internal transition, Eq. (22). The plane 
wave exponential would then be written as a product of an exponential term dependent on the CM coordinate, and 
the exponential term dependent on the relative coordinate, in complete analogy to Eq. (291. The explicit calculation 
of the transition amplitudes regarding the internal motion can be found in Ref. | 20 j . The main advantage of this 
treatment is that the effects of the electron spin are taken into account both directly and properly. Such a relativistic 
treatment does not change the scalar factors which couple atomic internal transitions to the electromagnetic field; they 
are proportional to Y(.- mt . An important modification is found in the polarization factor that in the non relativistic 
case is encoded in e TOph , and that now, distinguishes electric and magnetic effects through the presence of either 
electric or magnetic spherical polarization vectors [ 20 ]. 
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V. CENTER OF MASS TRANSITION INTEGRALS I CM . 

Since both the internal and center of mass states are simultaneously coupled to the electromagnetic field, controlling 
transitions between internal atomic states requires the proper regulation of center of mass transitions. During the 
last years, there has been impressive developments on laser cooling and trapping techniques, as well as an acquired 
capability to implement electric and magnetic fields on design. This has lead to the possibility of creating atomic 
samples in which each atom can be properly described by a wave function of precisely the form 


T(R,r reZ ) = y ^Cj 
ij 


(R)0 K rei (r rei ). 


(42) 


In particular, this applies to some experiments of few ions in Paul and Penning traps [2Tj and, very recently, to neutral 
atoms magnetically trapped in electromagnetic cavities |22| . In the case of ultra cold atomic gases in the degenerate 
regime, and for the condensate fraction, Eq. (421 can be used with $ k cm determined, e. g., by the Gross-Pitaevski 
equation for boson samples. 

In most other cases, environmental effects cannot be ignored and a description in terms of a density matrix for the 
center of mass state would be the most realistic option. The main reason behind this fact is the small value of the 
energy difference between the center of mass states of an isolated atom, compared to the interaction energy associated 
to collisions with other atoms and to the interaction with the background electromagnetic field. The implementation 
of semiclassical approaches where the center of mass motion is described in terms of Newton equations, while the 
internal dipole transitions are dealt within a quantum approach has probed to be adequate for the description of 
recoil effects in most standard situations [151 [2311251 . An alternative could be to take, as starting point, an idealized 
configuration in which each atom behaves as in empty free space with the center of mass stationary states taken as 
plane wave (PW) functions. Then the center of mass integral 1cm, Eq. (27), can be evaluated 


T FW ( k-lA 

J CMl K : K 0 ,^F 


,CM. .CM 


) = J(k — Ak cm), 


Ak, 


CM 


— L CM *CM 
— K F ~ K 0 • 


(43) 


The corresponding transition amplitudes, Eq. (261 (which is expected to be dominant), Eq. (30) and Eq. (361 would 
incorporate elementary recoil effects. The resulting transition rates should then be introduced in Bloch equations 
where other properties like collision and temperature effects could also be taken into account. Such a calculation is 
beyond the scope of this work. Nevertheless, in the following section, the consequences of using the center of mass 
factors i R the evaluation of the transition amplitudes is discussed for structured modes. 

Since ultra cold atoms are confined into a given region of space, in many cases, a reasonable description of the 
center of mass states is given by harmonic oscillator states centered at a position R 0 . Then, the center of mass wave 


function has an structure 4>* CM (R— Ro) and the transition integral Iqm, Eq. (27), takes the form 

k q r 


Icm( k; k 


CM. CM \ 
0 ! h F ) 


J d 3 R^lcM(R-no)e lk R ^ K cM 


(R Ro) 


_ gik-Ro 


d 3 R$*c M (R)e ik - ii $ K cM(R) 


(44) 


The exponential factor e ikH ° can also be written in terms of the wave vector k and the position of the center of the 


trap R 0 using again Eq. (371. Notice however that for the center of mass coordinate, ujRq/c may be greater than one. 


For the particular case of an atom confined in a potential that has either spherical symmetry or cylindrical symmetry 
around the z-axis the integral in Eq.(44l can be directly evaluated (see Appendix C). In the first case the dependence 


on k is through the spherical harmonics: 


t ST 

l CM 


(k; k\ 


CM. 

0 5 


Kp M ) — IcM {tC 2 /Mt^cmC 2 K\ 


CM. 


Kp M \£)Y ec 


.(k). 


Finally, if the center of was motion is described in terms of free spherical waves: 

!cm { k ; Kq M ; k f M ) = k$ m ; Hp M k). 


(45) 


(46) 


with A dB the de Broglie wavelength of the center of mass. The vector integrals I cm = ^klcM, Eq. (311, in both 
cases can be directly written in terms of vectorial spherical harmonics (Appendix B). 











VI. TRANSITION AMPLITUDES T 0 f FOR A CENTER OF MASS DESCRIPTION IN TERMS OF 

PLANE WAVES. 
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- \/^f - (Afc^)2)J+ i( Ak CM ;^;«^) 


If the center of mass factor Icm is given by Eq. (431, then 

TfS? = J dujS{cv)e~ iAt 5(Ak^ M 

For electric dipole (ED) transitions, 


iW*mAk™,Ak? M ) (47) 


T, 


,PW:ED 

OF 


1 


m p h 


J dojS(cj)e- iAt 6(Ak^ M - (Afc£ M ) 2 ) 


' (W cM e )v {A) ( l , l 1 , 0,1; to, to', 0, m ph )e* mph ■ e^¥^{Ak^ M , Ak° M ), (48) 

which is directly proportional to the angular spectrum as could be expected. This equation is a quantum analog in 
momentum space to the standard dipole interaction term d • E(R) used in most semiclassical descriptions of ED 
recoil effects. From this equation we observe that recoil effects are determined directly by the polarization and angular 
spectrum of the mode under consideration. This result has direct physical consequences, e. g., a Bessel mode [5] or a 
Laguerre-Gauss mode 110 can induce center of mass rotations around the main propagation axis of the beam, while 
a Weber mode has probed to be adequate for inducing a center of mass motion with parabolic symmetry on cold 
atomic samples [26]. 

For electric quadrupole (EQ) and magnetic dipole (MD) transitions the transition amplitude has the structure 


T 


PW:A 

OF 


= f duS(uj)e- iAt 6{Ak° M - sjuji/c 2 - (A/fcfM) 2 ) 


m ph ,m i 


• ^i]L uK , K '(k^/cM e )v [ f (lj\l,l]m,m\m 1 ,m ph )el lph ■ e^h\ mi (Ak M )¥^(Ak^ M , Ak^ M ), (49) 

where the EQ or MD character of the transition is encoded in the factor v^ which depends on both the mode 
polarization coupling parameter m p h and the effective angular momentum factor mi. In this case, the effective 
interaction term yielding recoil effects depends on the polarization and products of the angular spectrum of the mode 
with first order terms of k. In configuration space, these products correspond to first order derivatives of the electric 
and magnetic fields. 

Due to its frequent use in actual experiments, we discuss in more detail axial symmetric scalar Gaussian modes of 
waist wj_, 


1 Fw^(kj_) —W Wj _e 


—wi k , 


For scalar circularly polarized modes, the exact expression of the normalization integral is evaluated in Appendix A, 


^ 2 7T 3 


1 


|Wwx| 2 = fe*> -^[^ 2 +^ 


-4a 2 - 1 
8 a 3 


e 2a Erfi(v^2a)] 


i -i 


a = 


WWj_ 


(50) 


with Erfi the imaginary error function Erfi(,z) = eri{iz)/i. The factor fkv is a natural scale of the energy of the mode 
that could be used for the quantization of the EM field. For our study, it just provides a reference useful for comparing 
the effects of different structured modes with the same EM energy on atomic transiti ons rates. For Gaussian modes, 

different choices of the waist w_l modify the relevance of the Yi mi (Ak ) factor, Eq.(401; that is, for paraxial modes 


Yi,o(Ak CM ) 


1 


-\/47r ’ 


Y h± 1 (Ak CM ) 


0 


so that the internal transition rates which would be dominant are those with v ^ (1, l', 1,1; to, to 7 , 0, m p h) 7 ^ 0. In the 
antiparaxial regime k± ~ w/c, that is wj_ 0 , 


F 1>0 (Ak CM ) 


Y h± 1 (Ak CM ) 


e ±icpk C M 

\Z8n 


and the internal transition rates which would be dominant are those with v^\l, l', 1,1; to, m', ±1, m p h) 7 ^ 0. Notice 
however, that for a Gaussian mode the transverse momenta distribution is always centered around kx = 0 discarding 
them as the optimal choice for achieving the antiparaxial regime. A better option corresponds to using modes with 
a kx distribution that becomes maximum at a k± value different from zero, like the propagation invariant modes 
mentioned at Section II. 
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VII. TRANSITION AMPLITUDES T 0 f FOR ATOMS TRAPPED IN A SPHERICALLY SYMMETRIC 

POTENTIAL. 


If the atom is trapped in a spherically symmetric external potential of natural frequency (jJcmi the structure of 
the scalar center of mass term Icm hr the transition amp litude is given by Eq. (104), while that of vectorial internal 
transition factor J+ e[ has the structure given by Eq. (38). Since the product of two spherical harmonics can always 
be written as a superposition of spherical harmonics, ail the transition amplitudes Tof can then be written in terms 
of the basic integrals 


n 


( b-m p h',K ,) 
£,rri£ 


= / d 2 k x e. 


TTlph 


3 W#W(k ± )Y,, ro ,(k). 


(51) 


Similar arguments apply to the vectorial factor I cm and the scalar factors J^ el , so that these basic integrals also 
determine the transition amplitudes associated to the interaction Hamiltonians H 12 , Eq. (30) and H 13 , Eq. (36). 


In general, any integral H\ 


( b;m ph -K) 


is relevant for the evaluation of more than one transition amplitude. For instance, 


internal electric dipole transitions induced by Hu involve Yoo(k) which may be coupled to a center of mass transition 
with no rotational effects, that is Iq M = £p M and itiq M = m^ M . In such a case, the center of mass term involves 
a term Yoo(k), and the integral %o,o should be evaluated. A very different situation corresponds to an electric 
quadrupole transition induced also by Hu which involves internal Yi |TO (k) terms and a center of mass transition with 
a similar contribution (this transition could be accompanied by a change in the rotational wave numbers of the center 
of mass motion). The coupling of the two spherical harmonics will again have a contribution of the integral 'Ho.oj but 
with different internal and center of mass factors. 

In the following subsections we elaborate on the dependence of the on the electromagnetic parameters b 

and k for different EM modes. 


A. Spherical EM modes. 

Electromagnetic fields with spherical symmetry can be written in terms of vector spherical harmonics as described 
in Appendix B. If recoil effects are ignored in a radiative transition, the evaluation of internal transition atomic prob¬ 
abilities (including both electric and magnetic effects important in the relativistic realm) can be naturally described 
in terms of these modes ED]. In this section, contrary to standard calculations, recoil effects have been incorporated 
in the calculation. The integrals necessary to study the role for a given EM spherical mode on the 

coupled angular momentum changes in the internal and center of mass motion (as encoded in a with m p h 

= 0, ± 1) are 


n 


( b;ni ph ;jm) 


Vhuj 

27r 


d 2 k ± e 


TTlph 


■ Yg(k )Y tmi (k) = J dn k (rfl (k)) ^Y lmi (k). (52) 


Note that, the vector integrals I cm = Vk Icm, Eq. (311, are orthogonal to magnetic vector modes, the terms of the 
Hamiltonian Hj 2 that involve I cm would yield center of mass transitions just in the case that electric spherical modes 
are used to induce them. 

Using equations that can be found in Appendix B, it can be shown that 

v/ ^ [ r r n 

,j+l°me,m+m p h V d 


laJ (E;m ph ;jm), 

H e,m e 


I n 


( M;m ph ;jm ). 


27TC - 

+ j —1 time ,m-\-m p h VJ + 1 


3 + 1 1 3 

m + m ph -m ph -m 

3- 1 1 3 

m + m p h —m p h -m 

j 1 3 

m + m ph -m ph -m 


(53) 


here, E and M denote electric and magnetic spherical vector waves. The 8 factors indicate which spherical waves 
induce the coupled (internal and center of mass) atomic transitions. Taking into account the behavior of the 3j- 
symbols, in general, 



( b;m ph ;jm) 


< 


1 

2tt’ 


b = E,M. 


(54) 
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This result will be useful for the clear comparison between the integrals f or different selections of the 

electromagnetic modes. 

Though spherical EM modes give a natural scenario for the description of spontaneous internal transition rates, 
induced atomic transitions using these modes would require an efficient procedure for their selected generation; this is a 
non trivial task. Perhaps the simplest implementation could be achieved by the controlled excitation of predetermined 
EM modes within cavities with spherical symmetry. 


B. Gaussian modes. 

For axial symmetric scalar Gaussian modes of waist w_l, the integrals of the EM field that go along the polarization 
vectors e± necessary to evaluate atomic transition amplitudes are 


2 f i 

K ( tr£ ) = 8m,oMw ± y/(2t + 1)tt^ / dtPiW 1 - t 2 )e-“ V t, 

c Jo 


(55) 


while the integrals of the EM field that go along the polarization vector e z are 


a/(*>«) _ 

rL £,m e 


2 pi 

8 m , oA/w ± vW+ 1K^ / dfft(\/l-f 2 )e-“ 2f Vl-< 2 . 

c Jo 


(56) 


Some specific expressions are: 


K 


(±,wl) _ j\f /-U 2 1 e 
0,0 _ - /v wiV 7 ‘ 


M = o’ Wx) =-A/w x \/37r 


,2 

—uj 2 2 a — ^/7re~ aZ ErR(a) 


2 2 

uf,r x) =v wl ^ 2 a 2 ~ 3 + fT (3+a2) . 


L 2,0 


4a 4 


(57) 


These functions are plotted in Figure 1. It can be observed that it is necessary to consider modes with values of 
a = ww_l/c ~ 1 to observe a substantial increment of the transition rates with respect to the result obtained for 
paraxial beams. In general, the maxima of ; are achieved for waists wi smaller than the wavelength. As 

it could be expected, the functions decay as w_l increases faster than the functions 1 ^ do. In fact,for 

a » 1 


-i/(±,wjl) \/2(2l + 1) (2m)! prp 1 

^ ,0 7r 2 £ m\m'\ wi’ 


(58) 


m = m! = 1/2 if £ is even, and m = (£—l)/2 = m! — 1 if £ is odd. As a consequence for a scalar Gaussian mode, taking 
into account the normalization factor given in Eq. (50), all transition rates depend on the energy per unit transverse 
area of the EM mode in the paraxial limit. 

We have also evaluated the integrals T~L^ b ^ lph ’ v ' ± ' > for transverse electric and transverse magnetic Gaussian modes. 
The integrands differ from those of scalar beams by factors k±c/u ;, k±k z c 2 /u 2 and k\c 2 /uj 2 as a consequence of 




Eq. (hOk and Eq. (511. They also decay as w^ 1 as a » 1, but they yield maxima with higher values of 77) 
the scalar modes. In fact some of these maxima are beyond the spherical modes limit, Eq. (54). 


(K ’ Wi) than 


C. Propagation invariant modes. 

In Appendix A, the normalization of propagation invariant modes is worked out in several schemes. Two of them 
refer to the case where ideal propagation modes are considered. The third is the closest to the experimental realizations, 
it corresponds to work with a Gaussian superposition of ideal modes around a given k± : 


kj_ ,#c m :wj_ 




(59) 
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Figure 1: (Color online) (a)Illustrative basic integrals ’Hj'rn for scalar Gauss modes; and (b) illustrative basic 
integrals 'Hj'i, for vector TE and TM Gauss modes. They are plotted as a function of their waist w_l measured in 

terms of the inverse wave number a = cow± jc. 


with w, 1 k±. The relevant integrals for scalar propagation invariant modes can then be written as 


n l,mi L t,mi 


l(b, kj_) __ 

a e,m t ~ 


( 21 + 1 ) (t \me\)\ u- f 1 dt j^(b) p\ m e\ /./ 1 _ t 2 \ -CL 2 

2 (£+| TO ,|)l C 2 J 0 IVl t )e 


t = Kj_c/u 


1 


r 2TT 


qM = ^=I ^Me imev 

V27t Jo 


(60) 


The set of terms /()''_ is common to all propagation invariant modes, while (-){,)” depends on the details of the angular 
spectra t hat characterizes each kind of propagation invariant modes. 


2 a 


illustrates the absolute value of the integrals \hJ^ nj | = hj ^ for scalar modes as a function of the 


Figure 

dominant transverse wave number k± measured in units of u/c. In Fig 2b the vectorial analogues of these integrals 


are shown; their evaluation incorporate the terms e* h ■ e^ for b = TE, TM. Though both scalar and vector integrals 


depend on the waist of the Gaussian w ^ 1 used in Eq. (|59|, the same reported general structure is observed whenever 


w ± x k±. Note that there is a substantial increment of the hj^ n( values as the parameter k± deviates from the 
paraxial approximation. These occurs both for scalar and vector beams. 

The integrals 0^”^ have a very simple expression for Bessel beams: 


r\(B es sel:m) 
y ^'me 


1 


r 2n 


= l “ ne 


irrupt imetpk 


= 8- r 


(61) 


which is a manifestation of conservation of total angular momentum during the transition. Notice that this angular 
momentum can be originated from polarization terms like those in Eq. (10 ) or from orbital terms related to the 


angular spectra Eq. (191. They can also either be transferred to the center of mass or to the internal state of the 
atom depending the origin of the e lrnev± factor in the integral 0 m,. 

For Mathieu beams Q^ athleu - m '> can directly related to the A™ and B™ Mathieu coefficients in the standard 
notation|13j. They achieve a maximum value for q = 0 for m = me as could be expected since in this limit the elliptic 


coordinates reduce to circular cylinder coordinates, and Mathieu modes Eq. (20 1 become identical to Bessel modes, 
Eq. (191. For moderate values of q and as it increases, the off diagonal terms to 7 ^ me also increase. 


This means that Mathieu modes can be used to create superpositions of internal and center of mass states whenever 
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Figure 2: (Color online) Basic k± integrals, for (a) scalar propagation invariant modes (b) transverse electric 
and transverse magnetic propagation invariant modes as a function of the transverse wave number k± measured in 
units of w/c. The waist of the Gaussian factor was taken as lOc/w. 
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Figure 3: (Color online) Absolute value of the angular spectrum factor | = 0 mm / for Mathieu modes as 

a function of the continuous Mathieu parameter q = fn _l/2 with / the focal distance that defines the elliptic 

coordinates, and k± the transverse wave number. 


the corresponding internal, J re i and center of mass, Icm , vector or scalar factors have similar values for the mi values 
of interest. 

We now report the results for the absolute value of &i^' e ebe ' as a function of the continuous parameter a that 
defines the aperture of a Weber mode. In this case Q^ eber ' a ^ is much less than one, in general, with the exception 
of odd Weber beams with a nearby well identified discrete values. As a consequence, if a transition involving even 
Weber modes takes place, the resulting atomic states will be a superposition of many orbital angular momentum 
states. In general, due to the intrinsic central field basic structure of the internal atomic states, these means that the 
recoil effects of the transition on the atom will be better described not as simple rotational. This is consistent with 
the fact that Weber modes are not eigenfunctions of the orbital angular momentum but carry, as a natural mechanical 
property, the product of the angular momenta along the z axis and the linear momentum along a transverse direction 
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Figure 4: (Color online) 
modes, and (c) 


Absolute value of the angular spectrum factor \Q^ eber ' a ' > \ = Q^ eber ' a ' 1 (a) for even Weber 
for odd Weber modes as a function of the continuous parameter a. For all modes, 

0 (Weber:a) (Weber:a) 

me ~ 'O-mp 


p y , L~p y \ 27]. For internal dipole transitions induced by Hji, the center of mass factor Iqm will have a maximum 
overlap with the angular spectrum of a Weber beam if 

jj,OL In | tan ^k/2| 

Icm ~ 7 . . = j~- 

2y 7r| sin <^ k | 

This could happen if the initial and final wave functions of the center of mass are solutions of the Schrodinger equation 
of a free particle with boundary conditions with parabolic symmetry. One important question is whether or not the 
corresponding transition can be described in terms of a conservation law involving the atomic and the electromagnetic 
dynamical property L z p y . A first experimental study with that scope was reported in Ref. [26] . 


VIII. ATOMS IN A TRAP WITH CYLINDRICAL SYMMETRY INTERACTING WITH EITHER 

LAGUERRE-GAUSSIAN OR BESSEL BEAMS. 


For atoms in a cylindrical trap with the symmetry axis parallel to the main propagation axis of the light beam, 
according to the formalism developed in the first sections, the relevant integrals have the structure, 


C 


(b-K.) 


■' i , m i , m p h ; A m r »A*, <2 p , & 

For the Tq_f matrix element, 


= / d 2 k ± lZ Am ^(k ± a p )Z n °’<(a z J(uyc2 - k\) 


e m ph ' e 


(6) 


^ 6) (k±)lWk). (62) 


T 0 F = Vrel^ J dbjS(uj)d ernern ^(^/cM e -, K o ■ 


m p h ; Am s ,/i,a p ,a z 


(63) 


the parameters m p h, £ and mf are directly connected to the internal atom transitions as described in Eq. (381, that 


is m p h keeps the information about the vectorial coupling of the electron current to the EM field, £ = 0 for internal 
dipole transitions, and £ =1, m^=0,±l for electric quadrupole and magnetic dipole transitions. Higher multipole 
internal atomic transitions can be evaluated taking into account that, in terms of k± and <ys kj _, 


iwk) = N^ me e im ^VpWl - (fcic/u;) 2 ) 


(64) 


with the associated Legendre functions. The parameters Am#, /i, a p ,and a p in Eq. (621, encode the information 
of the center of mass transition matrix elements, Eq. (97||101 ) in Appendix C. The functions 1Z and Z exhibit an 
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exponential factor 


nz ~ e -“X/ 4 c 2 e -fcl(ai-a?)/4_ 


(65) 


For a spherical symmetric trap the equivalent exponential factor is recovered, Eq. (1031 Appendix C. The electromag¬ 
netic field structure is determined by the polarization vector and the angular spectrum (kj_). 

For Bessel modes the quantum numbers are the polarization b, the transverse wave number n± and the topological 
charge The evaluation of the integrals is direct, for scalar beams 


r (±;m B ,K±) _ • rsc [t r _ __r sr ] ^i,me 

{ -'m ph ,£,m e ;Am R ,/J,,a p ,a z ~ • /v kj_ ,ui [ u m ph ,±l"m ( ,-m B n. - f— - 1~^2 Um ph.,O u rn e ,m B ±l \ ^ £ 

y 1 — {K±C/ UJ) 


T,,0Om e ,m B ±l\ m R ,)i,a p ,a, 


.(«-l) 


( 66 ) 


while for vector beams 


p(TE;m B ,n±) — J\f’ 

'~'£,m£,AmR,ti,Oip,aL z J 


,U1 y/2 cui [ ^ rn P h ’ 1 ^ rrli ’ 1 - mB 3m ph ,-l8-m e ,m B + l] ^ Am R ,n,ot p ,a. z ( K -L) (67) 

K±Cy/l — K\c 2 /iiJ 2 


p(TM;m B ,K±) _ \fTM 

£,m£;AmR,fi,a.p,a z ^ k±,uj 


\ph 


Wl 


\ ^m ph ,l^m e ,l—m B + dm ph ,-l<5~m e ,m B + l] 

.2 , ,2 


K 2- U} S 

n u me ,—r 


c 


£,mi 

Amn,n,otp,<y : 


,(«l) 


with 


(«-0 = 2t TN e , me K Am ^( K± a p )Z^<(a z y/(u; 2 /c 2 - 4)7^ (^1 - («ic/ U )2). 


( 68 ) 


(69) 


These analytical expressions allow an optimization procedure of the parameters to achieve enhancements or inhibitions 
of the transitions rates. This is directly observed from the delta factors which determine the values of the Bessel beam 
topological charge for a given beam polarization that yield specific changes in the atomic quantum number; but it 
also applies for the selection of the kj_ parameter as illustrated below. 

Laguerre- Gaussian (LG) modes, usually studied in the paraxial regime fU 0 have an angular spectra given by 
Eq. (17). The vortex factor (k±wo/c) rn e zmcpk ^ leads to similarities and important differences between the integrals 


related to Laguerre Gaussian and to Bessel modes: while both share the same value of the integration on the tpk ± 
variable, LG beams become zero as k± —> 0 in a continuous way while for Bessel beams the variable k± takes abruptly 
the nonzero value n±. The angular similarity has as a consequence that both Bessel and LG beams can be used 
to enhance or inhibit transitions involving selected changes in the atomic (CM and internal) angular momentum 
as observed in Eq. ( 66 || 68 1 , and mentioned in the last section for CM motion with spherical symmetry. A detailed 
analysis on such a direction is reported, neglecting retardation effects, for Bessel beams in Ref. ||>] and for Laguerre 
Gaussian beams in Ref. (3]. In the latter case, the interaction electric Hamiltonian of the PZW formalism is used, free 
motion is assumed for the center of mass in the z direction, and the particular case p = 0 of the paraxial Laguerre 
Gaussian modes are studied. A quite interesting addition theorem for the factorization of the center or mass and the 
internal motion adequate for cylindrical symmetry is also introduced in Ref. . The allowed atomic and center of 
mass transitions for a given polarization and topological charge of the beams as described by Eqs.( 66p8| is consistent 
with the reported results in both references when the appropriate limits are taken. 

Let us illustrate an analytical procedure that can be used to optimize the selection of the w 0 (LG beams) and 
k _l (Bessel beams) parameters to enhance a given transition rate. LG transitions are determined by the elementary 
factors, 


e -^al/Ac 2 I" e -* 2 wl 

Jo 

The function e~ x w lg/ 4 ‘(-\/1 — x 2 ) s x r becomes zero at x = 0, 1, if s, r yf 0, and the factor (y/l — x 2 ) s x r achieves its 
maximum value at x 2 m = r/(r + s) < 1. If r = 0 and s / 0 (or s = 0 and r / 0) the maximum is achieved at x m = 0 
(or x m = 1). As a consequence 


/ 4 (\/l — x 2 ) s x r xdx, w 2 lg = ( a 2 — a 2 + w 2 )u> 2 /• 


(70) 


e -x 2 wla/i^i- x 2y x r xdx < (1 _ x^y^x 7 ^ 2 / e 


i, (l — x 2 W 2 x r ^ 2 1 — e - ™^/ 4 

l^xdx = l 1 Xm 1 e _ 

2 ° " 


<c/* 


and 


2 , , 2 r 1 2 2 , /- (l - T 2 W 2 -W 2 p-a ?^ 2 / 4c2 _ p-(a 2 +<" 2 )^ 2 / 4 c 2 

e-“<~ /4c / e- xw ^/\Vr^) s x r xdx< {1 Xm) Xm - 2 2 //i 2 - ( 71 ) 

Jo 2 ( a 2 - aj + w 2 )u) 2 / 4c 2 
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The maximum value of the latter term is one, and it is achieved for an elongated trap with a 2 = a 2 + w 2 —» 0, 

exhibiting and exponential decay on any direction of the ( a z , yj ot 2 p + w 2 ) plane, with independence on the values of 

the positive parameters r and s. As a consequence, transition rates are enhanced if the Laguerre Gauss beam is highly 
focused ( w 0 lj/c as small possible) and the trap is an elongated one with a 2 = a 2 + w 2 as small as possible. This 
is consistent with earlier observations |5] that predict that atomic transitions can be enhanced by locating the atom 
close to a light vortex. Note that for a p Cl, working with paraxial LG beams yields 


0 —utol z /Ac‘ 


r 1 2 2 ; ,_ Cl _ r 2 W2 r r /2 —cdV/4c 2 

/ e~ x w ^/\y/x - x 2y x r xdx < l 1 x ™> Xm e w 2 co 2 /4c 2 » 1 

J 0 2 w z oj z /Ac 2 


with the last factor being much smaller than one. 

For Bessel beams, and given values of the parameters a~ jP , the external parameter k±c/oj can be chosen to achieve 
the maximum value of the modulus of C. Calculations can be performed in an analytical way for the elementary term 


e -^/4cV*>y4 ( y^| )X) Xo = K±c/i0 , w 2 B = (a 2 p -a 2 z )co 2 /c 2 
which is also zero at x a = 0 and x a = 1 and achieves an extreme value for x Q given by 


(72) 


k b c 2 w B /2 + s + r =F y/{w 2 B /2 + s + r) 2 — 2 w 2 B r 
U} 2 w B 


(73) 


the sign is selected according to a 2 > a 2 or a 2 < a 2 . If a 2 = a 2 , the maximum is achieved at k]_c 2 /uj 2 = r/(r + s). 
These expressions give the optimal values of n± as a function of the exponents r and s, and the geometry of the trap 
in terms of the parameter w B . 


IX. CONCLUSIONS. 


In order to design experiments towards the control of atomic transition rates via laser light shaping, it is necessary 
to perform quantitative equitable calculations that allow the identification of the optimal parameters to inhibit or 
enhance such rates with respect to standard conditions where, in most cases, Gaussian laser beams in the paraxial 
regime are used. 

In this work we have shown that these calculations can be easily performed whenever a prior identification of the 
internal and center of mass states of interest is done; this follows from the fact that the structure of radiative atomic 
transition amplitudes can always be written in terms of: 


(a) An internal atomic transition factor. It has a vectorial character for the dominant interaction Hamiltonian Hu 
that results from the coupling of the electron current density to the electromagnetic field, and a scalar character 
for other terms in the non relativistic interaction Hamiltonian including one term in the coupling of the nucleus 
current density to the electromagnetic field, and the terms quadratic in the electromagnetic field. The vectorial 
factor Jr,; can be evaluated using different approaches to the description of the atomic dynamical space, for 
instance, a non relativistic approach as that implicit in Eq. (24 ) or a fully relativistic approach including the 
electron spin as the one given in Ref. |2flj . The central field structure of atoms makes that the dependence of 
the internal atomic transition factors on the electromagnetic wave vector k can be naturally expressed in terms 
of spherical harmonics Y^ m (k) in all cases. 


(b) An atomic center of mass transition factor that describes elementary recoil effects. Under standard conditions, 
the velocities P/Mt involved in the CM motion are much smaller than those associated to the internal motion 
Prei /Ad that reflects on smaller values of the matrix elements of the interaction Hamiltonian that couples the 
electromagnetic field to the nucleus current density than those related to the EM field-electron current density 
coupling. The atomic center of mass transition factor has a scalar character for the dominant interaction 
Hamiltonian Hji and for the terms quadratic in the electromagnetic field, and a vectorial character for the other 
terms in the minimal coupling interaction Hamiltonian. The vectorial factor I cm w & s shown to correspond to 
the gradient in wave vector space of its scalar analogue Icm- The center of mass transition factor can play 
an active role in the control of internal atomic transition rates. This has already been recognized in Ref. [5], 
where it was shown that the probability that the internal state of an atom acquires orbital angular momentum 
from a Bessel mode is maximum for an atom located close to the symmetry axis of the beam. For a center 
of mass Gaussian packet, the relevant parameter was found to be the ratio of the spread of the wave packet 
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to t he t ransversal wavelength of the photon. In the present formalism this correspon ds to taking Ro = 0 in 
Eq. (44) and tuu 2 /MtIjJcmc 2 ~ 1 to maximize the Icm{Tuj 2 /Mtucmc 2 ) value in Eq. (104). In the case of an 


idealized CM description in terms of plane waves the calculation of the elementary recoil effects is direct. The 
CM transition factors are necessary to perform, via their incorporation in the corresponding Bloch equations, 
realistic simulations of experiments at room temperature and moderate densities. 

(c) A factor that involves the EM radiation structure which is coupled to the internal transitions and recoil terms. 
The EM radiation is characterized in this term by its angular spectra ^f^u(k x ,k y ) and its polarization which, 
in general, also depends on k. 

The electromagnetic structure is completely encoded in the factor (c) that can be optimized to enhance or inhibit 
the transitions of interest. The general procedure has been exemplified for atoms with center of mass available states 
described by (i) plane waves, (ii) wave functions exhibiting spherical symmetry and (iii) wave functions for t rap s 
with cylindrical symmetry. In the second case, EM shaping is determined by the integrals given by Eq. (511; 

a quantitative comparative analysis between the transition rates for different EM fields was shown to be direct, as 
illustrated for spherical waves, Gaussian beams with different waists and propagation invariant beams with diverse 
angular spectra. 

Spherical modes provide an ideal scenario for the description of spontaneous atomic transitions between internal 
states. Their usage for the control of induced transitions is conditioned by the capability of creating spherical modes 
with predetermined characteristics. An upper bound for the important functions for spherical waves was 


found, Eq. (541. This value is also an upper bound for most of the corresponding integrals obtained for Gaussian 
(Fig. 1) and propagation invariant modes (Fig. 2). An interesting exception corresponds to transverse magnetic 
Gaussian beams inducing it transitions (e p /j = eo = e z ) in the antiparaxial regime, Fig. 2b. 

This study also allowed the recognition of the interactive role of polarization and EM angular spectra. For electric 
quadrupole transitions, induced by Hu, the selection rule A J = 0, ±1, ±2 allows changes of 2 h in the atomic angular 
momentum. This yields an implicit requirement of light modes with orbital angular momentum. These modes can 
be spherical (as worked out in Section VILA) or with other symmetries like Bessel modes or Laguerre Gauss modes. 
The plane waves decomposition in terms of spherical modes, Eq. (371, makes explicit their orbital angular momentum 
content. 

By changing the ratio between the transverse and longitudinal components in the dominant wave vectors yielding 
a structured beam, the relevance of predetermined internal states involved in a forbidden transitions can be modified. 
This is a direct physical consequence of the behavior of the functions TL < l , K ' r ^ M ' > that encode the polarization and angular 
spectra effects on transition rates. By evaluating these functions, it was shown that, in the paraxial limit and for 
scalar Gaussian beams, any multipole rate depends on the EM energy per unit transverse area. We also noticed that 
measurable increments of the transition rates for Gaussian beams with respect to this asymptotic value requires waists 
of the order of the wavelength. This is not a trivial task and, since Gaussian modes (as well as Laguerre-Gauss and 
Hermite-Gauss beams) have a transverse wave vector distribution centered around kj_ = 0, they are not necessarily 
the optimal choice for achieving that enhancement. 

For propagation invariant modes, we have noticed that they could be used to enhance transition rates in more 
favorable conditions since their approximate implementation involves a narrow wave vector distribution centered 
around kj^ ^ 0. 

The relevance of the angular spectrum on recoil effects was studied in general and illustrated for Mathieu, Bessel and 
Weber beams. Using a plane wave description of the center of mass motion, we showed that the EM angular spectrum 
can be directly imprinted in the atomic recoil effects for dipole transitions. In particular, our calculations confirm 
that, via dipole transitions, photons with orbital angular momentum can transfer this dynamical variable to the center 
of mass [Mi. For electric quadrupole and magnetic dipole transitions, the angular spectrum and polarization have 
diverse and interesting consequences on internal and center of mass transition rates. For those transitions, the effective 
interaction term describing recoil effects depends directly on the gradient of the EM field, so that semiclassical forces 
could depend on second derivatives of the electric and magnetic fields. It was also shown that Mathieu beams can be 
used to entangle atomic internal and center of mass degrees of freedom when the available states of the latter have 
spherical symmetry. 

Another important observation is that the formalism developed in this work can be adapted to other atomic 
systems such as molecules or nanostructures. In such cases, the particle that describes the transition of interest is not 
necessarily an electron. One could be interested, for instance, on the rotation or vibration of atomic nuclei. In any 
case, the symmetry of the internal available states should be taken into account for the most efficient calculation of 
the elementary integral containing the EM factor e* k r . 

If the atoms of interest are embedded in a material media, the electromagnetic field can be described using an 
equation similar to Eq. ([l]) by replacing the proper dispersion relation (which may depend on the polarization). 







19 


Evanescent waves would also require a modification on the Fourier factor e* k r to incorporate an exponential decay in 
the adequate direction. Then, the general formalism presented in this work would be properly adapted to describe 
transition rates in those systems. 
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Appendix A. Normalization of the EM modes 

The quantization of the electromagnetic field in terms of modes is usually based on the orthonormality requirement 

^ J d 3 r (E»* (r) • E<*> (r) + B £>* (r) • B$ (r)) = tw8{w - a j')S K , K >, (74) 

with the integration performed over the whole space SR 3 . For modes with polarization TE , TM this condition has a 
simple expression in terms of an integral of the modulus of the angular spectra, 

2tt 2 [ dV—^^(kL)*«(kx) = hu& K<K ,. (75) 

J LO OJ 


That is not the case for structured scalar modes due to the fact that the polarization vectors G^ 1 ) and are not 
orthogonal. This kind of modes are frequently used in experimental realizations where circularly polarized Gaussian 
beams impinge atomic samples. For scalar localized modes, Eq. (751 can be taken just as a normalization condition, 


i_ / ( B (±)' . E (±) + BW* ■ BW) = A f d 2 k _[l + W 


2 


k 2 c 2 . 


I'MiOl 2 - 


(76) 


In both Eqs. ( 75]|76 l, the integral in the left side must be performed within the circle of radius u/c in the kx space. 
For axial symmetric scalar Gaussian modes of waist wj_, 


=Mv 


,-wj k z . 


and circular polarization, the exact expression of the normalization integral valid beyond the paraxial approximation 
is 


1 

87 r 


EW* . E^ + BW* . B<±>) = \K 


2 ^ + ' /5¥ 


4a 2 - 1 
8 a 3 


e" 2 a 2 Erfi(v / 2a) 


= hio, 


(77) 


with a = wj_u;/c and Erfi the imaginary error function. In the case of Hermite-Gaussian modes of waists w XtV that 
guarantee the validity of the paraxial approximation, \k z \ ~ ui/c the polarization vectors and A 2 ' 1 are quasi 
orthogonal, and the integral over the circle of radius uj/c can be replaced to that over the whole space. Then the 
orthonormalization condition is approximately achieved by taking J\f = A WiWyfiw/2 m+n 7r 3 m!n!. 

Propagation invariant beams are non localized waves (the electric fields do not become zero fast enough as |x| —> 
oo)so that the orthonormalization requirement Eq. (741 gives rise to delta functions on the modulus of k±. Explicit 
results for TE and TM modes have been reported for Bessel[S], Mathieu|12| and Weber[T5j modes. 

An alternative for non localized waves, corresponds to evaluate normalization factors within a finite space region 


eW*.eW+bW*.B®) =hu. (78) 

For plane waves this normalization usually considers a parallelepiped of volume V = L 1 L 2 L 3 . This method has 
the advantage of being more intuitive and closer to experimental realizations. Taking this into account, we outline 
the main steps of this procedure for other propagation invariant beams. In that case, the natural space region for 
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performing the normalization condition is a cylinder with symmetry axis along the main propagation axis, radius 
R and length L. The normalization factor should then be written in terms of the angular spectra given above. A 
relevant result is: 


[ d 3 re i(k ' _k) r - 2ttR 2 L 

■M|kj_ — kj_ i?)i 

"sin {k' z — k z )L/2' 

'Vr,l 

k'j_ — k_L i? J 

1 (k' z -k z )L/2 


(79) 


Then, given a function of the form 


'Mr) = j d 2 k_L^ K (kj_)e zk r , k z = - k\, 


(80) 


iv R , 


d J r$;,(r)$ K (r)=2rfi / d 2 k^kx^k'JMkx) 


M|kj_ - k±|i?)i rsin(fc(, - k z )L /2 


kj_ — kj_|i? JL (k' z -k z )L /2 


(81) 


For L » c/uj , the last term is non negligible just for k' z ~ k z and, if we consider two modes with the same frequency, 
necessarily |kj_| ~ k' ± |. For propagation invariant TE or TM beams, the orthonormalization integral becomes 


1 

87T 


EW* ■ E<W + B?* • Bi b)N ) = 


.2 „2 


/ <i w [ <WJ(*v)/foW ,( ' /5 '‘ x ' R ' /1 ~ cos(w “ w)) 


47TW 2 


\/2k_l-R\/1 - cos(<£k' ^ M 


(82) 


For i? >> k^ 1 , 


od-Mv^kl-Rx/I - cos(w - M) v x , ft r\ - 7 -^ 

2R - -= - -► o(v2 k_l\/ 1 - cos(yv - v?k))- 

V2kj_Aa/ 1 - cos((/ 5 k ' - M 


(83) 


So that, 
1 


— jf (eW*.eW+bW*-bW) =nRL^ J dy* J dtp k f*,(<Pv)fMS((V2K X ^ ~ cos(^ - M) 


(84) 


For Bessel beams and Mathieu beams the integral over </j k is simply 2ir6 m i tm . For Weber beams, the parameter a is 
continuous and 


27r A(ol — a) In | tan c^k/2| ^ r oo 


/ rf< ^ k a I ■ 

/o 4-7 t| sin ipk 


— _ j e %(a-a‘)u _ 


2n 


= S(a - a'). 


(85) 


An third alternative valid for propagation invariant modes corresponds to defining them in terms of a Gaussian 
superposition of plane waves centered at kj_ instead of the idealized definition given by Eq. (181, 


ky) = k x = ^Jkl + kl ^ = atan (k y /k x ), 


( 86 ) 


with wj_ k±. In the limit wj_/r_l —> oo with finite n± one recovers the idealized definition. The elementary 
integrals necessary to evaluate the normalization factor for Bessel, Mathieu and Weber Gaussian beams 

are those described in the previous paragraph besides the integral 


dk±k±e 


—2v/ 2 L (k± — n,±) 2 _ 


2w l J W±K± 


dte 2t = 


2w_|_ AC_|_ 

4w^_ 


(87) 
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Appendix B. Spherical EM waves A complete basis for transverse electromagnetic fields in free space can be 
constructed in terms of spherical vectors. In the notation established by Eq. ([8]) , the corresponding modes are defined 
by (see, e. g., EQ]) 


<£?l( ka) ky-,w) = ^^ Yg(k), b = E,M, k = £i. 


k 

Iki 


In these equations, the superscript specifies the electric (E) and magnetic ( B ) spherical modes, and 


Y S-m(®k^k) = k X 


{E)( 


with 


„ * d 1 d 

V k = 0k WW" + ^k “—~7 - 


L <90k sin 0k ^k ’ 

Yj m (9 k> V’k) are th e spherical harmonics, and 

0k = cos 0k cos (p^ei + cos 0k sin ip^e 2 — sin 0k<33 

and 

<^k = — sinyj k ei +cos<pke2- 

It can be shown m that, 


Y Z = E ( Y i-)A«A, b = E,M 


A=0,±1 


( _ 1) i+m + A +1(Y g )A= _ CA ^( 


3 +1 1 3 \ y 

m + X -A -m 1 3+hm+x 


( 88 ) 

(89) 

(90) 

(91) 

(92) 

(93) 

(94) 


+ CA V^ I (m + A -A —inj ^“ 1 > m +A 

( _ 1)J+m+ A+i (Y (M) )A = _ CA% /2jTl ( m J + A E) (95) 

with Cq = i, c+ = —* and c_ = *. From these expressions an explicit relation for the basic integrals \'H\ E r ' Y:>n ‘ > | defined 


in Eq. (521 can be found, 
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Appendix C. Transition center of mass integrals for trapped atoms. 
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For the particular case of a confining potential that has cylindrical symmetry around the Z -axis we write 

1 VW- e- p2/2a2 pe- z2/2a Hp/a p ) lmcMl L^ CM \p 2 /a 2 p )H n ^z/a z )e xrncMipR ( 97 ) 


n p ,mcM ,n 2 


(R) = 


Jal(p+\mcM\)'. 


with a 2jP = ^'Kh/MTU> ZtP i and p = (n p + \mcM\)/2- The relevant integrals are 


p2n 

/ dp R e ik± p cos VR e iArn,RV>R = 2ni AmR J AmR (k ± p) 

Jo 
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JO 
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2 J~ V 

and M the confluent hypergeometric function. In particular 
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( 101 ) 


For the particular case of a confining potential that has spherical symmetry and natural frequency ujcmj the above 
results can be applied taking a = a z = a p = ttH/Mtujcm- Alternatively, we can write 


^ncM,loM,moM (R) — Rn C M,l (R)Y lm (e,<f>), 


( 102 ) 


where 


Rn^j(R) = N„ OMl{ R/V2 a )‘e-#/^Lii:\RP**), N naM , = ^ 

In that case it results convenient to use Eq. ( |37| , so that 
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The radial integral is 
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and the angular integrals can be written in terms of 3-j symbols yielding the coefficients c(l,m,l' 

Ic T M (h*j 2 /M t ujcmc 2 -, Kq M ; = c(l, m, l', m'; £ C m, tticm^cm^ 2 / m t^cmc 2 \n' CM li) 
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